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Introduction 


1.0  Introduction 

The  one  way  shape  memory  effect  (Buehier  and  Wiley,  1965,  McNichols  and  Cory,  1987, 
Schetky  and  Wu,  1991)  is  associated  with  twinning  induced  permanent  strain  recovery 
upon  heating.  Embedding  of  shape  memory  afioy  (SMA)  fibers  in  slender  flexible  bodies 
results  in  shape  changes  of  the  host  medium,  whenever  shape  recovery  of  the  SMA  takes 
place.  There  are  two  methods  of  embedding  SMA  fibers  in  flexible  structures  (Chaudhry 
and  Rogers,  1991).  The  prestrained  SMA  fibers  are  either  inserted  in  a  sleeve  and  then 
clamped  at  both  ends,  or  they  are  continuously  bonded  with  the  flexible  structure.  In  the 
first  case  the  result  is  a  point  follower  force  applied  at  the  ends,  while  in  the  second  case  a 
distributed  actuation  force  is  applied  from  the  fiber  to  the  flexible  structure.  This  paper 
deals  with  modeling  of  the  deformed  shapes  of  flexible  rods  with  continuously  embedded 
line  actuators. 

The  shape  change  of  a  cylindrical  rod  with  a  single  off-axis  embedded  SMA  fiber  has  been 
modelled  in  a  recent  paper  by  Lagoudas  and  Tadjbakhsh  (1992).  The  distributed  axial 
compressive  force  and  bending  moment  due  to  phase  transformation  in  the  SMA  fiber  was 
evaluated  using  an  approximate  shear-lag  model  and  the  deformed  shape  of  the  rod  was 
found  by  numerically  solving  the  equations  of  equilibrium  for  the  rod.  In  this  paper  we 
present  the  generalization  of  the  theory  to  a  flexible  rod  with  a  line  actuator  embedded  in 
an  arbitrary  location  with  respect  to  the  centroidai  axis  of  the  rod.  The  equations  of 
equilibrium  of  the  rod  and  actuator  are  derived  in  a  systematic  way  for  rods  undergoing 
large  3-D  deformations  (bending,  twist  and  extension).  An  example  is  presented  for  the 
case  of  a  SMA  fiber  actuator  undergoing  repeated  phase  transformations  between  the 
martensitic  and  austenitic  phases,  for  which  the  shape  memory  loss  and  the  reduction  in 
the  actuation  force  are  evaluated. 

2.0  Kinematics  of  Flexible  Rods  with  Embedded  Line 
Actuators 

2.1  Refarencs  Configuration  of  Rod  and  Lino  Actuator 

A  rod  is  described  in  the  reference  configuration  by  the  centroidai  line,  C,  that  connects  the 
centroids  of  its  cross  sections  and  by  the  reference  body  orthonormal  tnad  £.,  /  =  i,  z  3, 
where  £,  is  the  unit  tangent  vector  to  C  and  £,,  are  aligned  with  the  principal  axes  of  the 
cross-section  (Fig.  1).  The  actuator  is  geometrically  characterized  by  its  center  line,  i.e.  the 
tangent  vector  to  it.  We  will  assume  that  the  line  actuator  has  no  bending  resistance  and 
as  a  result  there  is  no  need  to  select  a  body  frame  dictated  by  the  principal  axes  of  its 
cross-sections.  As  it  will  become  obvious  later  on,  a  natural  choice  for  a  body  frame  for  the 
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actuator  is  the  orthonormal  triad  formed  by  the  tangent  vector  to  its  center  line,  A,  the 
principal  normal  and  the  binormal  to  A  (trihedral  basis). 

The  position  vector  x  of  the  line  actuator  in  the  reference  configuration,  with  respect  to  a 
fixed  Cartesian  coordinate  system,  is  given  in  terms  of  the  position  vector  x  of  the 
centroidal  line  C  of  the  rod  by 


X  =  X  +  P  (E0 1) 

where  p  =  dj:^  =  d^Ei  +  d^E^  is  the  position  vector  of  the  actuator  relative  to  the  centroidal 
line  C,  as  shown  in  Fig.  1 .  If  /  =  i,  2, 3  is  the  orthonormal  triad  of  the  fixed  coordinate 
system,  the  reference  body  triad  is  expressed  in  terms  of  the  direction  cosines  Ly  by 
E,  =  L/S)nj,  and  eq.  (1)  reduces  to 

X  =  (X,.(5)  +  d^(S)L.^(S))  n.  (EQ2) 


where  S  denotes  the  arc  length  along  C.  The  relationship  between  the  arc  length  Z  along 
the  actuator  A  and  arc  length  S  along  C  is  given  by 


1  =  ^  S 

dZ' dZ 


-2 


dp  dp  dp 
‘■'^^^y'ds^dsds^ 


(E03) 


The  above  equation  reduces  to 


dZ 


dd. 


dd^ 


i  =  J  (1  -  ^1^2  +  ^2^1 ) "  +  -  W  +  ^d,K,) 


dS 


^dS 


where  the  curvatures,  K,-,  i=1 ,2,3,  of  the  reference  shape  of  the  rod  are  defined  by 

dE. 
dS  ~ 


(E0  4) 


(EOS) 
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FIGURE  1  Reference  Configuration  of  Rod  and  Line  Actuator 


2:2  Dafomtod  Configuration  of  Rod  and  Una  Actuator 

Assuming  that  in  the  deformed  configuration  the  relative  position  of  the  line  actuator  witii  respect 
to  the  body  frame  remains  the  same,  the  position  vector  of  the  deformed  shape  of  the  actuator  is 
given  by 

X  -  X  + 4  (EQ6) 

X  -  {xfS)  +  dJ,S)l^^{S))n.  (EOT) 

where  x  is  the  position  vector  of  the  deformed  centroidal  line  of  the  rod,  c,  and  the  relative 
position  of  the  actuator  is  given  by  d  =  (Fig.  2).  The  deformed  body  triad,  1  =  1, 2, 3,  is 
expressed  in  terms  of  direction  cosines  Ijj  by  »  /y,<5)n^ . 


KlnMiMlic*  of  Floxlbl*  Rods  with  Emboddad  Lino  Actuator* 


FIGURE  2  Deformed  Configuration  of  Rod  and  Line  Actuator 


The  deformed  arc  length  z  along  the  actuator  is  found  in  a  procedure  similar  to  the  one  used  for 
Z,  with  the  only  difference  that  the  elongation  of  the  centroidal  line  c  has  to  be  taken  into 
consideration,  i.e. 


1  -  ^  ^  -  (4L\ 

~  dz  'dz  "  W 


((l  +  e)2  +  2(l+<r)f,a  +  ^.fe 


dS  dS' 


where  the  strain  e  is  defined  by 


dx 

55 


(1  +6)^3 


(Ea9) 


Eq.  (8)  finally  reduces  to 


Equations  ol  Equilibrium  for  noxIMo  Rods  with  EmbadcM  Una  Actuators 


dz 

dS 


dd. 


ddn 


(i-he-d^k2  +  d2k^)^+ -^2^3)  +  (^  +^1^3) 


dS 


with  the  curvatures,  kj,  1=1 ,2,3,  of  the  deformed  shape  of  the  rod  defined  by 

dS  ~ 


3.0  Equations  of  Equilibrium  for  Flexible  Rods  with 
Embedded  Line  Actuators 

3.1  Equilibrium  of  Lino  Actuator 

The  vector  form  of  the  equations  of  equilibrium  for  the  line  actuator  is 


The  internal  force  f  in  the  actuator  is  tangential  to  its  center  line,  i.e. 


(EQ  10) 


(EO  11) 


(EQ12) 


—  F^t  (BO  13) 

where  the  tangent  unit  vector,  r,  has  the  evaluation 


dz  ^  ddj  dd2 

( (  ^"dS  ~  ~  1  ^"dS  ^2^1  ^ 


The  distributed  force,  f ,  is  applied  from  the  rod  to  the  actuator  and  is  measured  per  unit 
undeformed  length  of  ttie  actuator.  If  the  trihedral  basis  is  used  to  resolve  the  vector  form 
of  the  equations  of  equilibrium  into  components,  the  following  set  of  two  equations  obtains 
(equilibrium  in  the  osculating  plane) 


dZ 


(EQ15) 


F^k+j^  =  0 


(EQ  16) 


The  actuator  curvature,  k,  per  unit  undeformed  length  in  eq.  (16)  is  defined  by 


Equations  of  EquUlbrlum  for  Floxibto  Rods  with  Embsddad  Lins  Actuators 


dz  '  dZ 

The  distributed  force  has  in  general  the  following  resolution  in  the  osculating  plane 

where  the  principal  normal  vector  is  defined  by  (Sokolnikoff  and  Redheffer,  1966) 

=  i^  - 

~~kdZ  kdS^dS^ 

3^  Equilibrium  of  Rod 

The  vector  form  of  the  equations  of  equilibrium  of  a  flexible  rod  acted  upon  by  a  distributed 
force,  /,  acting  at  a  distance  d  from  the  centoid,  is  (Love,  1944,  Tadjbakhsh  and  Lagoudas, 
1992)' 

dF  r 

^+f,  =  0  (E<1») 

If  we  resolve  the  forces  and  moments  in  the  body  coordinate  system,  we  have  the 
following  set  of  six  equations 


(EQ17) 


(EQ  18) 


(EQ  19) 


dF. 

dM, 

35  +  «.;/ -<!+«>  ^,73^;  +  Vt  =  0  <“  “> 

The  explicit  form  of  the  equations  of  equilibrium  is  also  given  here  for  later  reference 

dPx 

^  +  A:2F3  - /:3F2 +/j  =  0  (eq24) 

dF^ 

TS  +*3f,-*,f3+/2  =  0 


(EQ25) 


Equations  of  EquUlbriutn  for  Raxlbia  Rods  vwith  Embsddsd  Uns  Actuators 


dF. 

55  “^2^1  ■‘■/b  “  ® 

dM. 

^  -  ( 1  +  c)  +  <^2/3  “  ^ 

dM^ 

^  ‘•rk-^M^  -^jA/3+  (1  +e)F^  “^1/3  ~  ® 
dM^ 

+k^M2-  ^2^  1  +  ^1/2  ~  ~  ® 

In  the  above  equations  Ft(s),  F2(s)  and  F3(s)  are  the  two  shear  forces  and  the  axial  force 
in  the  rod,  respectively;  k^,  k2,  k3,  are  the  two  curvatures  and  the  twist  of  the  deformed 
shape  of  the  rod  per  unit  undeformed  length,  respectively;  K-i,  K2,  K3,  are  the  curvatures 
and  the  twist  of  the  undeformed  shape  of  the  rod  in  the  reference  configuration, 
respectively.  The  moments  of  inertia  of  a  cross  section  of  the  rod,  with  respect  to  the 
principal  axes,  are  denoted  by  li^  and  l22t  while  the  polar  moment  of  inertia  of  the  cross 
section  is  denoted  by  I33. 

The  distributed  force,  (,  applied  on  the  rod  from  the  actuator  is  connected  through 
Newton’s  third  law  to  ^e  distributed  force,  f ,  applied  from  the  rod  to  the  actuator  (Fig.  3)  by 

r  dZ  rn 

f  = 

Notice  that  to  be  able  to  find  the  components  of  /  in  terms  of  the  components  of  the 
distributed  force  ^  =  f,i  +f„n  must  be  resolved  in  the  body  coordinate  system. 
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FIGURE  3 


Distributed  Forces  Applied  to  the  Rod  and  Line  Actuator 


If  linear  constitutive  relations  are  used  the  three  bending  moments  and  the  axial  force  are 
related  to  the  curvatures  and  extension  by  the  foiiowing  equations  (Tadjbakhsh  and 
Lagoudas,  1992) 


Afj  =  £/j  j  (/Cj  -  (EQ31) 

M2  —  f/22  (^2  *”  ^2^  (EQ32) 

Afj  =  G/32  (ikj  -  ATj)  {EQ33) 

F3  =  AEe  +  E  {I  l22^2^^2~  ^1)^  (6034) 


E  and  G  are  both  material  constants  relating  the  difference  in  the  curvatures  with  the 
corresponding  moments  and  can  be  approximately  taken  equal  with  the  Young's  modulus 


Equation*  of  Equillbriutn  for  Raxibl*  Rod*  with  EmtMddad  Una  Actuator* 


and  shear  modulus  of  the  material  of  the  rod.  Eqs.  (24-29)  have  to  be  solved  for 
F^,  Fj,  e,  k^,  Atj.  *3.  If  the  rod  is  inextensible,  i.e.,  e=0,  eq.  (34)  is  not  valid  any  longer  and 
becomes  an  unknown  variable  replacing  e  in  eqs.  (24-29). 


The  connection  between  the  curvatures  and  the  body  angles  of  rotation  (p,ei,  (pje, 
(Kane,  Likins  and  Levinson,  1983),  which  bring  the  body  triad  from  the  fixed  basis  n,,  n^, 
to  the  deformed  configuration  e^,  e^,  is  explicitly  given  by  (Tadjbakhsh  and  Lagoudas, 
1992) 


k,  = 


= 


^  cos<P2COsq>3-H^  sin(p3 

Ts  cos(p2Sin(p3 

*3  =  S  +3S  ““*^2 


(EQ35) 


(EQ36) 


(EQ37) 


Similarly  the  connection  between  the  reference  cunratures  and  the  body  angles  of  rotation 
<!»,£,,  OjEj,  OjEj,  which  bring  the  body  triad  from  the  fixed  basis  n,  to  the  reference 
configuration  E^,  e,  is  explicitly  given  by 

COSO2COSO3 -I- ^  "^811103  (EQ38) 

*^2  =  cosOjSin^j  (MM, 


^3 

ds 


-H-j-  sinO., 
ds  2 


(EQ40) 


The  position  vector  of  the  centroidal  line  of  the  rod  is  given  in  both  the  deformed  and  the 
reference  configurations,  respectively,  by 

5  5 

=  J(l+^)  l^^dS  +  x^iO)  =  J(1  +  c)  sin(p2tfS-hXj(0)  (eq41) 

0  0 
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S  S 

X2  =  J  ( 1  +  e)  l22^S  +  ^2(0)  =  “  J  ( 1  +  ^)  sincp  j  cos(P2^iS  +  jc2(0) 
0  0 

5  S 

=  j  (\  +  e)  l^^dS  +  x^iO)  =  J(1 +^)  cos(pjcos(p2^i5  +  i3(0) 
0  0 

S  S 

=  jl^^dS  +  X^CO)  =  Jsin<I>2^/5  +  X,(0) 

0  0 


(EQ  42) 


(EG  43) 


(EQ44) 


S  S 

X2  =  jL32i/5  +  ^2(^)  “  -  JsinOjC0sO2^^  +  -^2(^) 

0  0 

5  5 

X3  =  Jl33JS  +  X3(0)  =  Jcos<I>jCOS<I>2^^  +  ^3(0) 

0  0 

In  the  case  of  an  inextensible  rod  (e=0),  eqs.  (41-43)  satisfy  the  constrain  ^  =  1  ■ 


4.0  2-D  Bending  of  a  Rod  by  an  Off-Axis  Line 
Actuator 

For  2-D  bending  about  we  assume  that  jf,  =  at,  =  0  9,  =  9,  =  0.  The  position  vector 

of  the  deformed  configuration  of  the  actuator  is  given  by 

X  =  X  +  d^{S)e^  (EG  47) 

and  the  tangent  vector  to  its  center  line  has  the  evaluation 

dd^ 


ddy^ 

J(l+£-rfl^2)  +  W) 


(E04«) 


2<D  Bending  of  a  Rod  by  an  Off*Axia  Una  Actuator 


If  we  define 


tan(|>  = 


dd^ 

Is 


1  +  e  -dyk2 


(EQ  49) 


the  tangent  and  principal  normal  vectors  of  the  line  actuator  take  the  form 

£  =  sin(|)ej  +  cosij)^^ 
n  =  cos(l)|e|  -  sin(|>f2 


(EQSO) 

(E051) 


The  above  formulation  leads  to  the  following  evaluation  of  distributed  forces  and  moments 
on  the  rod 


/,  =  -(cos(<|))yJ  +  sin(4.)y^)^ 

/j  =  -(cos(4i)yf-sin(it>);^)^ 

<”2  =  ^2jk‘‘/k  =  ‘'l  (1»>?  -  (f fs 

where  from  eq.  (4)  we  have 


(E052) 


(EOS3) 


(E0  54) 


^  = 
dS  1 


dd,^ 

2  .  “"1 


(E055) 


4.1  Lina  Actuator  Placed  Parallel  to  the  Centroidal  Line  of  a  Straight  Rod 
The  results  presented  in  the  previous  section  simplify  in  the  case  of  an  actuator  embedded 

dd, 

parallel  to  the  centroidal  line  of  an  initially  straight  rod.  Since  ^  =  o,  eqs.  (50,51)  simplify 
to 

£  =  ^3  (EQ56) 
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n  - 

(EQ57) 

and  the  curvature  of  the  actuator  from  eq.  (19)  becomes 

(EQ58) 

Finally,  eqs.  (52-54)  reduce  to  (note  that  (j)=0  and  d2/dS=1) 

V 

tl 

(EQ  59) 

> 

II 

(EQ60) 

mj  = 

(EQ  61) 

If  the  above  simplifications  are  taken  into  account,  the  equilibrium  equations  for  the 
actuator  (15, 16)  take  the  form 


II 

(EQ62) 

F-tj  =  /, 

(EQ63) 

while  the  equations  of  equilibrium  for  the  r<xl  (24  -  29)  reduce  to 

dF. 

(EQ64) 

dF^  ^pa 

(E065) 

35^(1 +«)F,-d,f^  =  0 


(EQ66) 


Eqs.  (64-66)  are  the  same  with  the  ones  reported  in  Lagoudas  and  Tadjbakhsh  (1992), 
except  for  the  term  k^F',  which  is  not  included  there.  This  term  becomes  important  for  large 
deflections  with  large  bending  curvature,  while  it  can  be  neglected  for  cases  with  small 
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resulting  bending  curvature.  The  kinematic  relationships  for  the  deformed  shape,  i.e.,  eqs. 
(35-37,  41-43),  likewise  reduce  to 


S 

=  |(. 

0 

S 

=  /(. 

0 

The  results  derived  in  this  section  will  be  used  to  model  the  response  of  a  rod  to  an 
embedded  SMA  fiber  upon  attempted  shape  recovery  of  the  SMA  fiber  due  to  the 
martensitic-austenitic  phase  transformation. 

5.0  2-D  Bending  of  a  Rod  Induced  by  an  Embedded 
Shape  Memory  Alloy  Fiber 

Let  a  SMA  fiber  be  prestrained  while  in  the  martensitic  phase  and  embedded  in  a  flexibie 
rod  and  bonded  with  the  rod  material  in  a  temperature  below  Mf  (Schetky  and  Wu,  1 991 , 
Pfaeffie  et  ai..  1992).  Upon  heating  above  Af,  the  martensite  will  transform  to  austenite  and 
the  initiai  ineiastic  strain  would  be  recovered,  were  it  not  for  the  constraint  provided  by  the 
rod.  For  a  flexible  rod  with  an  off-axis  embedded  SMA  fiber,  such  as  the  one  shown  in  Fig. 
4,  the  evaluation  of  the  actuation  force  based  on  a  shear-lag  model  has  been  given  by 
Lagoudas  and  Tadjbakhsh  (1992).  For  completeness  we  summarize  here  the  key  results. 


L  —  _ ^ 

^2  -  dS 


e)  sin(p2^5-i-ij(0) 


e)  costp^^*^  +  ^3(0) 


(EQ67) 


(EQ68) 


(EQ69) 
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The  shear-lag  equations  of  equilibrium  are 


2 

+— a,,  =0  0<S<L 

dS 


oa  3  1 

5;  Pa<'-<P-<'l 

dc  2 

jc - =0  O^S^L 


(EQ  70) 

(Eori) 

(EQ72) 


Eq.  (70)  gives  the  force  equilibrium  in  the  SMA  fiber  in  the  axial  direction,  eq.  (71 )  is  the 
force  equilibrium  in  the  shear  layer  of  thickness  (p-d^)  in  the  radial  direction  and  eq.  (72)  is 
the  force  equilibrium  in  the  surrounding  rod  material  in  the  axial  direction.  The  constitutive 
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equations  for  the  stresses  in  the  fiber,  shear  layer,  o,,,  and  rod  material,  o,  are 
respectively  given  by 


dw. 

—  r  (  “  _  c^ma\ 

-  ^57 


(EQ  73) 


(EG  74) 


(EG  75) 


where  Wa(S)  is  the  eixial  displacement  of  the  SMA  fiber  and  w(r,S)  is  the  axial  displacement 
of  the  surrounding  rod  material.  E  and  G  denote  the  moduli  of  the  rod  material,  while  is 
the  Young’s  modulus  of  the  fiber  in  the  austenitic  phase. 


Combining  eqs.  (70-75),  together  with  the  condition  of  overall  equilibrium  of  forces  in  the 
axial  direction,  the  final  equation  for  the  force  in  the  SMA  fiber,  r*  =  is  found  to  be 
(Lagoudas  and  Tadjbakhsh,  1992) 


where 


-  aV  = 


(EG  76) 


c-sma 


E 

-  pI 


a2  = 


2G 


(EG  77) 


(EG  76) 


We  will  refer  to  as  the  actuation  force,  because  it  is  the  driving  force  in  eqs.  (64-66)  and 
gives  rise  to  the  distributed  force  applied  to  the  rod  through  eqs.  (62, 63).  As  eq.  (73) 
suggests,  e'**  is  the  stress  free  recoverable  transformation  strain  when  the  SMA  fiber  fully 
transforms  into  the  austenitic  phase.  Imposing  boundary  conditions  F®=0  at  S=0  and  S=L, 
the  following  solution  obtains  for  the  actuation  force 
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F«(5)  =  1  - 


-aS  -a  {L-S) 
e  +  e 


I  +  e 


-aL 


(EG  79) 


If  the  above  evaluation  of  the  actuation  force  is  substituted  into  eqs.  (64-66),  these 
equations  can  be  solved  for  F^,e,k^,  and  the  deformed  shape  can  be  evaluated  from  eqs. 
(67-69).  We  have  solved  equations  (64-66)  numerically  using  the  shooting  method,  in  a 
similar  way  as  described  by  Lagoudas  and  Tadjbakhsh  (1992)  for  various  boundary 
conditions,  geometric  and  material  parameters. 


5.1  Simply  Supported  Rod  with  an  Off-Axis  Emboddsd  SMA  Fiber 

As  an  illustrative  example,  and  to  present  a  method  for  finding  the  deformed  shape  of  a  rod 
with  an  embedded  SMA  fiber  undergoing  repeated  phase  transformations,  we  report  the 
case  of  a  simply  supported  rod,  shown  in  Rg.  5,  with  an  off-axis  embedded  SMA  fiber. 


FIGURE  5  A  Simply  Supported  Rod  an  Embedded  SMA  Fiber 


The  boundary  conditions  for  the  simply  supported  case  are  given  by 


xj(0)  =  ^2(0)  =  ^3(0)  =  0 


(EQM) 
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Afj  (L)  =  ^2  “  ^3  ~  ^  (EQ«i) 

The  initial  stress  free  straight  configuration  is  assumed  to  occur  while  the  SMA  fiber  is  in 
the  prestrained  martensitic  phase  before  any  shape  recovery  takes  place.  When  heated 
above  the  phase  transformation  temperature,  Af,  the  SMA  fiber  transforms  into  the 
austenitic  phase,  attempting  to  recover  its  original  shape  before  prestraining.  The 
deformed  shape  of  the  rod  has  been  evaluated  for  material  and  geometric  parameters 
shown  in  Fig.  5  with  initial  prestrain  of  the  SMA  fiber  =  -4%.  We  have  assumed  that 
the  austenitic  phase  behaves  linearly  elastic  with  E/Ea=1/20.  The  stress  state  in  the  SMA 
fiber  is  indicated  schematically  as  state  1  in  Fig.  6,  and  the  resulting  deformed  shape  and 
curvature  are  plotted  in  Figs.  8  and  9,  respectively.  For  the  numerical  solution  of  eqs.  (64- 
66)  the  rod  has  been  assumed  centrally  inextensible,  which  is  a  valid  assumption  for  SMA 
fibers  with  small  diameter  compared  to  the  diameter  of  the  rod. 
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If  the  prestrain  of  the  SMA  fiber  is  large  enough,  the  actuation  force  may  reach  the  yield 
limit  of  the  austenitic  phase,  o^,  which  is  actually  the  starting  point  of  stress  induced 
martensitic  transformation  (Schetky  and  Wu,  1991).  The  condition  for  yielding  is  given  by 

=  "Pa°a 

The  l.h.s.  of  the  above  equation,  normalized  by  is  plotted  in  Fig.  7  below.  The 

normalized  value  on  the  r.h.s.  is  aporoximately  equal  to  the  ratio  j/  .  We  can 

therefore  conciude  for  this  case  from  Fig.  7  that,  in  order  not  to  have  yielding  at  any  point 
along  the  SMA  fiber,  the  yield  strain  of  the  martensite  must  be  at  least  0.6  of  the  prestrain, 
or  2.4%. 

FIGURE  7  Normalized  stress  in  ths  SMA  fiber 


SigmaSMA 


If  yielding  occurs  in  the  SMA  fiber,  we  assume  that  shear-lag  model  is  still  valid  and  the 
actuation  force  can  be  found  from  eq.  (79)  but  with  the  replacement  of  eq.  (77)  by 

C‘'  =  Vp«[£^ '“•» 

The  above  equation  implies  a  perfectly  plastic  model  for  the  austenitic  phase,  as  shown  in 
Fig.  6,  which  is  an  experimentally  verified  valid  assumption,  due  to  the  formation  of  the 
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compliant  stress  induced  martensite  (Sato  and  Tanaka,  1988).  The  above  procedure 
requires  an  iterative  solution  scheme  for  the  evaluation  of  the  bending  curvature  in  eq.  (83). 

After  the  rod  cools  down  below  Mf,  and  the  SMA  fiber  returns  to  its  martensitic  phase,  the 
rod  does  not  return  to  the  initial  straight  configuration,  but  it  has  some  residual  curvature 
due  to  residual  stress  development.  To  evaluate  the  actuation  force  in  this  case  we 
assume  that,  after  the  SMA  fiber  transforms  to  the  martensitic  phase  in  a  volume 
preserving  transformation,  it  is  stretched  by  the  deformed  rod  until  a  new  equilibrium  is 
established.  Since  the  martensitic  phase  is  much  more  compliant  and  with  a  lower  yield 
stress  than  the  austenitic  phase  (Schetky  and  Wu,  1991),  it  yields  giving  rise  to  new 
prestraining.  To  model  the  constitutive  behavior  of  the  martensite  we  use  an  elastic 
perfectly  plastic  model  as  shown  in  Fig.  6,  where  the  cool-down  process  is  indicated  by  the 
path  1-2.  The  assumption  of  a  perfectly  plastic  model  is  reasonable  because  the  inelastic 
strain  is  mainly  due  to  twinning  stretching  and  to  a  lesser  extent  to  dislocation  motion  with 
a  very  small  amount  of  hardening  (Schetky  and  Wu,  1991 ). 


The  actuation  force,  F'”,  in  the  martensitic  phase  is  found  by  the  shear-lag  model 
described  by  an  equation  similar  to  eq.  (76)  with  the  solution  given  by 


F'”(5)  = 


1  - 


I  +  e 


with  the  following  evaluation  for  the  constants 


2G 


E 


(p  -d^)/pj  P  -P^J 


^  m 


(EQM) 


(EQ85) 


(EQ86) 


where  Em  is  the  stiffness  and  Om^  is  the  yield  stress  of  the  martensitic  phase.  The 
deformed  shape  of  the  rod  corresponding  to  the  martensitic  phase  is  shown  as  curve  2  in 
Fig  8  and  the  bending  curvature  in  Rg.  9.  For  the  martensitic  phase  we  have  used  E/Em=l/ 
6andOm^/Em=0.01. 
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To  model  subsequent  cycles,  the  transformation  strain  has  to  be  corrected  to  include 
residual  stresses  in  the  SMA  fiber  after  transformation  from  austenite  to  martensite  at  the 
end  of  process  1-2.  The  new  transformation  strain,  is  given  by  the  initial 
transformation  strain  reduced  by  the  residual  elastic  strain  of  the  SMA  fiber,  i.e.. 


^sma  _  ^sma  _ 


f  F' 


m 


2r.  “^1^2 


'm 


Substitution  from  eq.  (84)  into  the  above  formula  results  in 


£•^'”^(5)  =  +  - j —  1- 


psma 
nt 


a  m 


I  +  e 


-pL 


Far  away  from  the  ends,  the  above  formula  reduces  to  a  simple  one 


(EQ87) 


(EQ88) 


r  (jy 

^sma^L^  =  ^  (tow 

nt 

utilizing  eq.  (88)  for  the  shape  memory  loss,  ttie  deformed  configuration  and  bending 
curvature  of  the  rod  is  evaluated  and  shown  as  cun/es  (3)  in  Figs.  8  and  9,  respectively. 
Subsequent  cycles  will  not  influence  the  residue  stress  state  and  the  rod  will  be  deforming 
between  states  (2)  and  (3). 

As  a  final  note  and  to  bring  some  connection  witti  tho  composite  rod  theory,  we  rewrite 
eqs.  (64-66)  in  the  following  form 

dF. 

-I- /c2  (F-j -I- F^)  =  0  (EQW) 

jL(Fj  +  f") -tjf,  =  0  (ta.i) 

+  (l+«)f,  =  0  team 

Away  from  the  ends  of  the  rod  and  for  vanishingly  small  diameter  of  tiie  SMA  fiber,  we  have 
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E 


= -ko^Ee 

a 


sma 


'a  P‘-P 


(£093) 


If  we  consider  ttie  actuation  force  as  an  internal  force  for  the  composite  rod-SMA  fiber 
system  due  to  the  phase  transformation  strain,  we  could  absorb  the  actuation  force  in  the 
constitutive  equations  (31-34),  i.e., 


M2  =  £/22  (*2  -  £2)  +  ‘'l  «o«) 

F3  =  AEe■vE(I\^k^(,k^-K^)  +I^k2{k2-K2))  ieqss) 

in  which  case  the  equations  of  equilibrium  become 

0  (EQ  96) 

0  (EQ  97) 

=  0  (EQ  98) 

For  small  curvature  and  extension,  equations  (94-98)  reduce  to  the  classical  composite 
beam  equations  (Allen  and  Haisler,  1985).  The  shear-lag  model,  therefore,  in  the  case  of 
line  actuators  is  compatible  with  the  composite  beam  theory  and  in  addition  takes  into 
account  the  end  effects. 


I 

dS  2^3 

dh 

5s  -*2^1  = 


dMo 


I 
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6.0  Inverse  Problems 

Equations  (64-66)  can  also  be  thought  of  as  differential  equations  for  the  actuation  force, 
F®,  provided  that  the  bending  curvature  is  known,  i.e.,  a  desirable  deformed  shape  is 
specified  a  priori.  In  the  case  of  inextensible  rods,  the  system  of  equations  (64-66) 
becomes  linear  in  F-j,  F3,  and  F^  and  it  can  easily  be  integrated.  Before  we  go  into  the 
solution  of  the  inverse  problem  in  general,  we  give  the  solution  to  a  special  case  below. 

If  we  eliminate  F3  between  eqs.  (64)  and  (65),  in  the  process  F^  gets  also  eliminated  and 
we  end  up  with 


dS{k2dS  j 


+  ^2^1 


=  0 


(EQ99) 


which  has  the  general  solution 

=  CjCOS(P2  + C2Sin(P2  (eqioo) 

where  Ci  and  C2  are  arbitrary  constants.  SubstitiJtion  of  F^  into  eq.  (64)  yields 

F3  =  -F^  +  CjSin<P2“^2‘'®^^2  (eqioi) 

To  obtain  M2  wo  substitute  oq.  (100)  into  oq.  (66)  which  reduces  to 
^^2  dF^ 

“^155  =  -(1 +e)  (Cj 00892  +  ^2 sin(P2)  (eqio2) 

Integration  of  the  above  equation  yields 

M2 =  -  CjiCj- C2ij  +  Cj  (EQ103) 

where  C3  is  a  constant  of  integration.  Considering  moment  equilibrium  of  a  segment  of  the 
rod  we  arrive  at  the  following  result 

M2(5)  =  d^F^iS)  -  P3  (jcj(L)  -X|(S))  +  (x^{L)  - x^iS))  +  T2  (eqio4) 

where  f  =  p^n^  +  is  the  externally  applied  force  at  S=L  and  r  =  is  the  applied 
moment  at  S=L.  If  sufficient  boundary  data  are  known  eq.  (104)  can  be  utilized  to  correlate 
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the  actuation  force  with  the  deformed  shape.  In  the  case  of  a  free-free  end,  Fi=0  at  S=0 
and  S=L.  Hence  from  eq.  (100)  Ci=C2=0.  Also  M2=0  and  F®  =0  at  S=0  and  S=L  and  from 
eq.  (103)  63=0.  Eqs.  (101)  and  (103)  simplify  in  this  case  to 

F3  =  -F^  ^2  = 


The  above  simple  solution  indicates  that  if  a  desired  shape  is  sought,  the  required 
actuation  force  is  proportional  to  the  difference  of  the  bending  curvatures  in  the  deformed 
and  reference  states. 


For  more  general  boundary  conditions,  the  following  procedure  can  be  followed  for  the 
solution  of  the  inverse  problem.  If  we  assume  e  -  0  (centrally  inextensible  rod), 
differentiating  eq.  (66)  twice  and  substituting  from  (64,  65)  leads  to 


A. 

dS 


[dS  ^^2 


2  ,  dm  > 

^^2^  '^dS^dS  ^^2^  +fc2"‘2  ~  ® 


(EQ  106) 


where  m2=-dif3.  If  we  substitute  for  M2  from  eq.  (32)  and  define,  ^ ,  the  above 
equation  yields  ” 

^ (^2  ■  ^2)  /^2)  ^2^  (^2  ■  ^2)  +  ^ ^/^2)  ^2^2  =  ^ 

This  equation  is  a  third  order  nonlinear  o.d.e  for  the  bending  curvature  ^2-  Notice  that  while 
eq.  (104)  was  derived  by  integrating  the  equations  of  equilibrium  (64-66),  eq.  (107)  is 
obtained  by  differentiating  them. 

For  the  solution  of  the  direct  problem,  where  the  distributed  actuation  force  is  given  and  the 
induced  deformation  is  sought,  eq.  (107)  or  equivalently  eqs.  (64-66)  are  solved 
numerically  using  the  shooting  method  (Lagoudas  and  Tadjbakhsh,  1992).  Eq.  (107)  can 
be  primarily  utilized  to  solve  in  closed  form  two  inverse  problems: 

1 .  Inverse  Problem  (a):  Find  m2  for  given  k2  and  K2  along  the  rod. 

2.  Inverse  Problem  (b):  Find  K2  for  given  k2  and  m2  along  the  rod. 

For  the  first  inverse  problems  eq.  (107)  is  a  linear  second  order  o.d.e.  for  m2  with  variable 
coefficients.  The  homogenous  solution  is  given  by 
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m2  =  Z>j  sin  +  Z)2COS 
^0  ^  ^0 


y 


(EQ  108) 


The  particular  solution  can  be  found  using  the  method  of  variation  of  parameters.  If  the 
initial  curvature  is  zero,  in  which  case  eq.  (107)  reduces  to 


.  dm2  .  (d}k2  1  ? 

4;(tc  =  -4c\^  /k^  +  ^k. 


dS^dS  "'2^  •  2"‘2  ^s[dS 
We  can  write  the  general  solution  of  eq.  (109)  in  the  following  compact  form 


(EQ  109) 


m2  =  -|j^(5)cos  (0(S)  -  0(5))  JS  + Dj  sm0(5)  +  D2Cos0(5) 


'0 


where  g(S)  and  6(8)  are  given  by 


(EQ  110) 


g{S)  = 


^  ^2  1,3 

. 

dS  "^2^2 

(EQ  111) 

(EQ  112) 

5 

=  \k2dS 

(EQ  113) 

0 

The  two  constants  Oi  and  D2  are  found  by  satisfying  appropriate  boundary  conditions.  As 
an  example,  for  a  simply  supported  rod.  •L<S<L,  eq.  (110)  reduces  to 


L  .  S 

^2  =  ^J«(S)co8(e(S)-e(L))JS-^(L)ji^g-J«(S)cos(0(«-0(S))d5 


(EQ  114) 


As  special  cases  we  report  certain  closed  form  solutions  for  the  first  inverse  problem.  If  the 
final  bending  cunrature  is  given  by 

5 

^2  ~ 


(EQ  115) 
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the  required  distributed  bending  moment  to  produce  the  above  curvature  from  the  initial 
straight  configuration  is  found  to  be 


-  *^0’'  •  s  . 

"■2=  2r""^"2Z> 

A  parabolic  curvature 


requires  a  linear  distributed  bending  moment  given  by 


(EQ  116) 


(EQ  117) 


(EQ  118) 


and  a  hyperbolic  curvature 


(EQ  119) 


requires  a  distributed  bending  moment  given  by 


XS  XS 
sech  (-^)  tanh  (-^) 


(EQ  120) 


in  the  above  equations  Kq  is  a  constant  measuring  the  amplitude  of  the  required  curvature. 
The  value  of  the  parameter  X  measures  the  spread  of  the  cun/ature,  which  a  pproaches  a 
delta  function  as  X  becomes  large. 


A  particular  case  of  the  first  inverse  problem  arises,  when  the  deformed  configuration  of 
the  rod  is  required  to  be  straight,  i.e.,  k2=0.  Eqs.  {64-66)  reduce  in  this  case  to 

.  ‘‘>‘2 

^2  =  ^  (E0121) 

The  above  equation  suggests  that  the  integral  of  the  distributed  bending  moment  is  the 
required  initial  curvature  to  produce  a  straight  deformed  configuration. 
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To  solve  the  second  inverse  problem  we  introduce  the  new  dependent  variable 
\  =  and  eq.  (107)  takes  the  form 


^  dA  j  dm2 

55^55 “^^2^  "^^2^  ~  ~  ^55^55  +^2^2^  (EQ122) 

The  above  equation  has  the  same  form  with  eq.  (109)  and  can  be  solved  in  exactly  the 
same  way  with  eq,  (109)  of  the  first  inverse  problem. 
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